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It is widely recognized that the main difficulty in designing devices which could process informa¬ 
tion using quantum states is due to the decoherence of local excitations about a ground state. A 
solution to this problem was suggested in pj, relying on (non-local) topological excitations, struc¬ 
turally protected against local noise. However, a practical implementation of this proposal using 
special Landau levels in fractional quantum Hall effect systems (FQHE) [2| has proven elusive, while 
accessible FQHE states are theoretically not optimal because their representations in the Hilbert 
space of states are not dense. We propose using a different physical system (cold Fermi atoms), 
whose semiclassical dynamics is described by a hyperelliptic function in the Sklyanin formalism. 

The homological structure of the complex curve corresponds to representations of the braid group, 
with the action of Hecke operators leading to singularities detectable in the semiclassical oscilla¬ 
tions. We argue that, for a fixed genus of the hyperelliptic curve, the Richardson-Gaudin pairing 
Hamiltonian problem is the singular limit m —» oo of level-fc sl( 2), with k + 2 = Sr ^ +1 —¥ 0, so that 
the level k is admissible in the sense of Kac and Kazhdan )3], but the corresponding Hecke algebra 
is a q— deformation of the symmetric group with fixed q = e"^ 4 , asm-)' oo. 


1. Introduction Two decades ago, the introduction of 
the first quantum algorithm |dj led to an intense activity 
aimed at finding physical systems whose quantum states 
could be used in a reliable way to process information. 
In a quantum system, information is encoded in quantum 
bits (or qubits), rather than a classical bit. Unlike a clas¬ 
sical bit, which can only hold a value of 0 or 1, a qubit 
can be 0, 1, or any irreducible superposition of both. A 
qubit is realized by way of a two state quantum system. 
A popular approach to development of a qubit has been 
to make use of Rabi oscilations. Rabi oscilations are an 
external induced sinusoidal time evolution of a quantum 
system between its two states, typically by way of op¬ 
tical or electromagnetic stimulation. This behavior has 
been fundamental to many approaches in realization of 
a working quantum computer. The Rabi oscillations can 
be geometrically represented through use of the Bloch 
sphere. The classical states of 0 and 1 correspond to 
the poles of the sphere. As the system evolves in time 
the state vector oscilates between the poles of the sphere 
by a variety of paths. However, in spite of experimental 
successes designing various types of qubits Q, a success¬ 
ful implementation of quantum algorithms is still elusive, 
mainly because of the fundamental effect of decoherence 
which characterizes all devices based on local quantum 
excitations. Decoherence is the irreversible loss of a quan¬ 
tum systems ability to process information; it arises when 
a quantum system becomes coupled to its environment 
(which is unavoidable if we wish to control the evolution 
of the system). This coupling results in the informa¬ 
tion previously contained in the system being lost to the 
environment rendering further computation impossible. 
Decoherence is represented on the Rabi sphere by the 
state vector rotating around the equator of the sphere, 
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in essence being neither 0 nor 1 nor a superposition. This 
effect, studied under various assumptions on the type of 
noise [d] describing the coupling of the system to an exter¬ 
nal bath, is very important because, for devices based on 
local excitations (or quantum oscillations associated with 
single-state quantum dynamics), noise and useful signals 
are indistinguishable. To address this important issue, 
A. A. Kitaev and collaborators jlj have proposed using 
non-local, topological excitations, which are basically un¬ 
affected by most types of noise from coupling with the 
environment. An essential ingredient in the construction 
of these special excitations is the fact that they require 
the simultaneous creation (or annihilation) of two exci¬ 
tations at distinct points zi,z 2 - As explained in jl], the 
main advantage of working with such states comes from 
their non-locality: the likelihood that noise would cre¬ 
ate two excitations simultaneously at Z \, z 2 is completely 
negligible, so that the system is protected against a large 
class of perturbations, namely local random fluctuations. 

2. Topological excitations from Quantum Hall states 
The specific physical model which seemed to allow a di¬ 
rect implementation of Kitaev’s proposal is the Fractional 
Quantum Hall Effect (FQHE), more precisely the par¬ 
tially filled Landau level at filling factor v = 5/2. Fol¬ 
lowing [/], we briefly recall the essence of this theory: 

(i) the ground state of the v = ^-filling fraction FQHE 
Landau level is given by the PfafHan construction 


^Gs(zi, ■ • ■ ,Z2n) ~ Pf 



II(‘* “ ^) 2 ’ 

k<l 


(i) 


(ignoring single-particle factors), while an anyon 
excitation (obtained by removing two particles, i.e. 
by creating two holes at 2 : 1 , 2 ) corresponds to 


dQ/iOi, 22 ; 23 , • • •, z 2n ) ~ Pf 


zi - z 2 

ZlZ 2 (Zi - Zj) 


n (zk-zif, 

k<l 
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up to single-particle factors, where Pf stands for the 
Pfaffian of an even-dimensional, skew-symmetric 
matrix, Pf(A) = ytDet(A); 

(ii) the two-hole excitation changes under a permuta¬ 
tion of its arguments as 

'& 2 h(z 2 ,zi;z 3 ,.. .,z 2n ) = e 2 ^ L ^2h(zi,z 2 ;z3,.. .,z 2 „), 

obeying the anyon statistics with k = 4 (k = 1, 2 
correspond to bosons and fermions, respectively). 

Unfortunately, another feature of this theory makes it 
theoretically not optimal for achieving quantum compu¬ 
tation, and prompts current research for finding similar 
non-commutative anyon theories: specifically, it belongs 
to a class of theories whose states form irreducible rep¬ 
resentations for the sl(2) k Wess-Zumino-Novikov-Witten 
models (they are those with level k = 1,2,4), which are 
however not dense in the Hilbert space of physical states 
[sl]; this constitutes a fundamental limitation. 

3. Alternative realizations of anyonic theories To ad¬ 
dress the problem, we make the important remark that 
the Pfaffian state described in the previous paragraph 
can also be realized in other systems, physically differ¬ 
ent from the one provided by FQHE. Therefore, it may 
be possible to realize the same non-local topological ex¬ 
citations, while circumventing the representation-theory 
complication. To begin, note that the effective quantum 
theory corresponding to anyonic excitations with k = 4 
is the sl( 2)fc =2 WZNW theory 

Fixing a level k £ N of the representation gives the 
central charge of the theory as c = 3fc/(fc+2), and leads to 
a number of possible types of physical excitations, whose 
commutation relation is of anyonic type, with fractional 
statistics. Moreover, the primary field g(z, z) satisfies the 
equations (Knizhnik-Zamolodchikov) 


[( k+2)d z -J a T a ]g(z,z ) = 0, J a {()g(z,z ) = -^—g(z,z), 

S z 

with r a the generators of the Lie algebra su( 2). 

It is relevant to note the correspondence between 
the Kac-Moody algebras sl( 2)* and the quantum group 
lA q {sl( 2)), where the quantization parameter q is a root 
of unity, q = exp[*7r/(fc + 2)]. This can be regarded as 
a manifestation of the Riemann-Hilbert correspondence 
(between monodromy of solutions of linear differential 
equations and the singular points of the equations). Al¬ 
gebraically, this correspondence identifies the centralizer 
of sl(2)k to the Hecke algebra H(q ), a q —deformation of 
the symmetric group, and providing representations for 
the braid group in the corresponding number of elements. 

f. Cold Fermi gases in the fast-pairing limit The 
physical system we propose to use instead of quantum 
Hall devices consists of cold fermionic atoms, trapped 
in an optical lattice, and driven to create Cooper pairs 
through the Feshbacli resonance CM3- It is well- 
known [18]-[24j that a proper model for the dynamics 


of this system is the Richardson-Gaudin pairing model 


31|, which we briefly review here. It describes a sys¬ 


tem of n fermions characterized by a set of independent 
one-particle states of energies ei,l £ A. Each state l has 
a total (spin) degeneracy di = 2, and the states within 
the subspace corresponding to l are further labeled by an 


internal quantum number, s =t, Let c\ s represent the 
fermionic creation operator for the state (Is). Using the 
Anderson pseudo-spin operators jUJ (quadratic pairing 
operators), satisfying the su(2 ) algebra 


[tl tf] = ±6ijtf , [t+, tj] = 2<v,r;. (2) 

the Richardson pairing Hamiltonian is given by 


Hp = J2 - 9 J2 2eit I ~ 9 t+ ■- (3) 

ig A i,v IgA 


where t = t/ is the total spin operator. It maps to 
the reduced BCS model 

H = e P 6 P,^P>^ - 9 Cp t cLp^c_kic kt (4) 

p,cr p.k 


by replacing the translational degrees of freedom by ro¬ 
tational ones, where I 6 A = ennumerates the 

one-particle orbital degrees of freedom, while s = f, 4- i n_ 
dicates the two internal spin states per orbital (di = 2). 
The pairing Hamiltonian can be decomposed into the lin¬ 
ear combination 


H P = 2 ^ eiRi + g 

IgA 



IgA 


(5) 


At a fixed value of the component t 3 of the total angu¬ 
lar momentum, the last term becomes a constant and is 
dropped from the Hamiltonian. The operators Ri (gen¬ 
eralized Gaudin magnets [3jJ ) are given by 


Ri = tf 



ti ■ t v 
ei - ev 


( 6 ) 


These operators solve the Richardson pairing Hamilto¬ 
nian because [20j they are independent, commute with 
each other, and span all the degrees of freedom of the sys¬ 
tem. Richardson showed [23,123 that the exact At—pair 
wavefunction of his Hamiltonian is given by applica- 

tion of operators b' k = Y^i 2il - ek to vacuum (zero pairs 
state). The unnormalized At— pair wavefunction reads 
'kfl(ei) = IlfcLi frfc|0). The eigenvalues e*, satisfy the self- 
consistent algebraic equations 

~ = z2 --• ( 7 ) 

9 ^ k ek ~ e P I 2ei - efc 

As an integrable model, the Richardson-Gaudin system 
was solved based on the methods of [33}, in [d ll i ni. 
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5. Richards on- Gaudin model as singular limit of Wess- 
Zumino- Witten theory In this section we recall the re¬ 
lation between the Richardson-Gaudin model and sin¬ 
gular 517(2) Chern-Simons and WZW models fd~2l. Rdl . 
These relations stem from the study of quantum states 
of the Chern-Simons theory with gauge group 517(2) 
on the manifold T 2 x I and in the presence of Wil¬ 
son lines {zi} x R, i = 1 [HJ, corresponding to 

the one-particle spectrum of the Richardson Hamilto¬ 
nian, Zi = 2 ei. The torus T 2 has modular parameter 
t = Ti + iT 2 , 71,2 G R. Thus, T 2 = C/Z + rZ. The 
quantum states of this theory are known to satisfy the 
Knizhnik-Zamolodchikov-Bernard (KZB) equations 

W C s = 0, V = (V Tl V Zi ), (8) 


6. Semiclassical solution of the Richards on-Gaudin 
model The solution starts from the Lax operator 


£(A) 


2 

g 


V3 + J2 

2=0 


Si ■ a 
A — €i 


a{ A) 6(A)' 

c(A) —a( A) ’ 


( 10 ) 


where a a , a = 1, 2, 3 are the Pauli matrices, and A is an 
additional complex variable, the spectral parameter. Let 
Uk, k = 1,..., n — 1 be the roots of the coefficient c(A). 
Poisson brackets for variables 5“ read 

{5“, 5f } = 2 e afh S%8 jk (11) 

The Lax operator m defines a Riemann surface (the 
spectral curve) T(y, A) of genus g = n — 1, through 


where V is the flat KZB connection, 

Vr = (k + 2)d T + H 0 ( Zi ), V 2i = ( k + 2)d Zi +H i {z i ), (9) 

and Ho, Hi are elliptic versions of the Gaudin Hamilto- 
nia (and depend on the torus parameter r). The pa¬ 
rameter k is related to the level of representation of the 
affine algebra sl(2)k- By a limiting procedure, the system 
degenerates into the Richardson-Gaudin (R-G) problem: 
by taking the limit r —»• ioo, the torus degenerates into 
a cylinder, and the elliptic Hamiltonia Hi, Ho reduce to 
a trigonometric limit. Upon rescaling such that the set 
{zi} is collapsed into the origin, the cylinder becomes the 
complex plane with punctures at {zi}, and the rescaled 
operators become the rational Gaudin magnets ([6)) . 

The R-G model is formally retrieved in the singular 
limit k + 2 —>• 0, which imposes a careful analysis of 
the semiclassical solution in order to identify the proper 
choice for the fractional levels k as c —>• — oo. This anal¬ 
ysis will be carried out in greater detail elsewhere, and 
we discuss here only the main conclusions. Choosing the 
sequence of fractional levels k = 2 to G N, or equiv¬ 
alently k + 2 = 8n ^ +1 , c = |(1 — 8 to), q = e l7r / 4 , the corre¬ 
sponding representation of the braid group H ( q ) remains 
fixed as to — y oo. The fractional levels are admissible in 
the sense of Kac and Kazhdan Q, and they are under¬ 
stood in the sense of jJzj]. This generalized notion of 
representation level settles in particular the issue of irre¬ 
ducible representations of the modular group SL( 2,Z). 
As noted in [47], the initial interpretation of fractional 
levels for sl( 2) [4£j led to a finite number of irreducible 
representations of 57(2, Z), which is clearly insufficient, 
as the set of automorphic forms on the Poincare upper 
half-plane H, invariant under 5L(2,Z), of integer weight 
in 2N is described by the union of the Berezin kernels 
from the tower of von Neumann algebras B(H n ), n > 1, 
where H n = L 2 (H, y n ~ 2 dxdy). In contrast to the origi¬ 
nal Kac-Wakimoto approach to fractional level represen¬ 
tations of sl( 2) (which leads to a finite number of irre¬ 
ducible representations of the modular group), the gen¬ 
eralization [49| provides an infinite number. 


V 
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Q(A) = det£(A) 



( 12 ) 


where P( A) = ~ e 0- 

The equations of motion become 


2iy{uj) 

Ylj^iiUi-Uj)’ 


ij~ = J~ 


n 

gj 3 + 2y ^e k -u 
k =1 


Here, u = -2 H u i) = °- 

From the equations of motion, it is clear that knowl¬ 
edge of the initial amplitude of and of the roots {iq} 
is enough to specify the n unit vectors {S,}, for a given 
set of constants of motion {Ri} given by the classical 
limit of Gaudin Hamiltonia. The Dubrovin equations 
of motion are solved by the inverse of the Abel-Jacobi 
map, as we explain in the following. We begin by noting 
that the polynomial Q{ A) has degree 2 n, and is posi¬ 
tively defined on the real A axis. Therefore, the curve 
T(?/,A) has n cuts between the pairs of complex roots 
[E 2 i -i, E 2i ],i = 1 , 2 ,...,n, perpendicular to the real 
A axis. The points ut belong to n — 1 of these cuts, 
Ui € [l? 2 i-i, E 2 i\,i = 1,..., n — 1. These g = n — 1 cuts 
allow to define a canonical homology basis of T, consisting 
of cycles {a,, /%}, i = 1,..., g. With respect to these cy¬ 
cles, a basis of normalized holomorphic differentials { 04 } 
can be defined, through 


M2 = A®" 1 —, Mij = 

y 


m, u> = m 1 fi. 


(13) 


The period matrix Bij = jp.oJi is symmetric and has 
positively defined imaginary part. The Riemann 6 func¬ 
tion is defined with the help of the period matrix as 


9(z\B) = e 2ni{ntz+ i ntBnS >. (14) 

n£zs 


The g vectors consisting of columns of B and the ba¬ 
sic vectors e*. define a lattice in C 9 . The Jacobian variety 
of the curve T, is then the g —dimensional torus defined 
as the quotien J(T) = C 9 /(Z® + BZ 9 ). The Abel-Jacobi 
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map associates to any point P on T, a point ( g — dimen¬ 
sional complex vector) on the Jacobian variety, through 
A(P) = ui. Considering now a g— dimensional com¬ 
plex vector of points {Pfc}, k = 1,..., g on T, defined up 
to a permutation, we can associate to it the point on the 
Jacobian 

g 

z = a(P)=J2 A ( P ^ + K ’ ( 15 ) 

fc=1 

where K is the Riemann characteristic vector for T. 

The map m suggests that we now have a way to 
describe the dynamics on T by following the image point 
on the Jacobian. Given a point on the g— dimensional 
Jacobian 2 = (£i,..., Cn-i)> we can find an unique set 
of points {Afc}, k = 1,... ,g on T, such that 2 = a( A), 


and 9(a(P) — z\B) = 0. The system evolves in time 
according to the point z{t) 

Ck = ic k , 1 < i < g - 1, Cn-1 = *(c n — 1 + t), (16) 

where {c*,} is a set of initial conditions, such that z 0 = 
z(t = 0) = a(c), and c is the set of initial conditions 
for positions of A on T. Together with the initial condi¬ 
tion which determines the initial amplitude of J~, this 
set completely determines the functions m(t), J~{t). The 
application of inter-twiners of the Hecke algebra H ( q ) will 
be reflected in the homological structure of the Sklyanin 
solution for the semiclassical R-G problem as the coales¬ 
cence of a pair of branch points (followed by their sub¬ 
sequent resolution), which can be translated directly at 
the level of the solution J~(t). 
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